Coexistence of spin- 1/2 and spin-1 Dirac-Weyl fermions in the edge-centered honeycomb lattice 
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We investigate the properties of an edge-centered honeycomb lattice, and show that this lattice features both 
spin- 1/2 and spin-1 Dirac-Weyl fermions at different filling fractions / (/ =l/5,4/5for spin 1/2 and / = 1/2 
for spin 1). This five -band system is the simplest lattice that can support simultaneously the two different 
paradigmatic Dirac-Weyl fermions with half-integer spin and integer spin. We demonstrate that these pseudo- 
relativistic structures, including a flat band at half-filling, can be deduced from the underlying Kagome sub- 
lattice. We further show that the signatures of the two kinds of relativistic fermions can be clearly revealed 
by several perturbations, such as a uniform magnetic field, a Haldane-type spin-orbit term, and charge density 
waves. We comment on the possibility to probe the similarities and differences between the two kinds of rel- 
ativistic fermions, or even to isolate them individually. We present a realistic scheme to realize such a system 
using cold atoms. 
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I. INTRODUCTION 

In the solid state, matter is typically organised into crystal 
structures. The mathematical models for describing different 
materials are consequently based on lattices where the elec- 
trons are trapped in periodic structures. The understanding 
of the equilibrium and transport mechanisms in such systems 
also forms our knowledge of many fundamental effects and in- 
deed technological applications of today. Recently emergent 
phenomena such as quasi-relativistic effects in non-relativistic 
settings have proven to be important in this respect. Most 
notably graphene 1 , together with topological insulators 2 and 
cold atoms in optical lattices 3-8 , are prominent examples of 
this. The emergent quasi-relativistic fermions in graphene and 
topological insulators are two-spinor massless Dirac fermions, 
but it has also been suggested that in more exotic lattices, such 
as the ^3 lattice 9 and the line-centered-square Lieb lattice 
10-12 , the emergent massless fermions are in fact pseudospin- 
1 objects which also involve a flat band. Also Dirac-Weyl 
fermions with arbitrarily large spin have been studied based 
on fermionic atoms trapped in optical superlattices 13 ' 14 . 

In this paper we explore a different direction: can we find 
a single setup where several different kinds of Dirac-Weyl 
fermions can coexist in the same lattice? This is an intrigu- 
ing question. Such a "material" should have remarkably ver- 
satile properties as far as density dependent effects are con- 
cerned, as we will show in this paper. Our initial efforts 
to solve this problem is motivated by the study of Lieb lat- 
tices 10-1 2 , where additional lattice sites on the edges of the 
square lattice give rise to a flat band. From this inspiration, 
we expect that when introducing additional lattice sites to the 
edges of the standard honeycomb lattice 1 , which we refer to 
as the edge-centered honeycomb (ECH) lattice in the follow- 
ing, a flat band should also emerge 15 , thus giving both spin- 
1/2 and spin-1 Dirac-Weyl fermions. This is indeed what we 
have found in this study. We also note that the next-nearest- 
neighbor (NNN) hopping in this ECH lattice produces the 
well-known Kagome lattice 16 , thus the ECH lattice interpo- 
lates between several well researched lattices, such as hon- 



eycomb, Lieb and Kagome lattices. In fact, by using an in- 
triguing mapping (see Appendix A), we find the band struc- 
ture of the ECH lattice is completely determined by its un- 
derlying honeycomb and Kagome sublattices, revealing the 
deep connection between the ECH, the honeycomb, and the 
Kagome lattices. Furthermore, we investigate the response 
of the system to perturbations such as a uniform magnetic 
field, a Haldane-type spin-orbit coupling 17 , and a charge den- 
sity wave 16 (CDW). We demonstrate how these perturbations 
allow us to probe the similarities and differences between the 
two kinds of relativistic fermions. This five-band model turns 
out to be a minimal model that can support simultaneously the 
two different paradigmatic Dirac-Weyl fermions at the lowest 
spin level where spin- 1/2 Dirac points coexist with a single 
spin-1 Dirac point crossed by a flat band. 



II. THE MODEL AND ENERGY SPECTRUM 

We are interested in the properties of a fermionic gas 
trapped in an ECH lattice whose unit cell contains 5 inequiv- 
alent sites, which are labeled by T = 1, . . . , 5, as illustrated in 
Fig. 1 (a). A physical realization of this system could be 
achieved by trapping fermionic atoms using six lasers, which 
divide the plane into six sectors of 60°. The corresponding 
configuration of the laser light fields can be chosen as 

Ei =£(0,l)e toai , E 2 =E(V3/2,l/2)e- ikx - a \ 

E 3 =E(V3/2,-/2)e ikx - a \ E 4 = E(0, -l)e~ ikx ^ , 

E 5 = E(-V3/2, - 1 /2)e ikx ^ , E 6 = E(-V3/2, 1 /2)e~ ikx ^ , 

(1) 

where ai = (1,0), a 2 = (-1/2, \/3/2), a 3 = (-1/2, -y/3/2) 
are the three nearest-neighbor (NN) vectors of the underly- 
ing honeycomb structure (i.e., the red sites in Fig. 1 (a)). 
In this case, the intensity profile I(x,y) = |E to t(-X,;y)| 2 from 
the total electric field E to t = E; produces a potential land- 
scape as shown in Fig. 1 (c). Alternatively, one can en- 
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visage using Spatial Light Modulators for shaping the inten- 
sity of a light beam such that the desired minima create the 
ECH lattice (see for instance Whyte and Courtial 18 and refer- 
ences therein). Also nano structured lattice potentials for two- 
dimensional electron gases 19 can be considered for creating 
the ECH lattice. 

For sufficiently deep lattice sites we can use the tight- 
binding approximation, which results in the non-interacting 
Hamiltonian 

ffo = *I>fo, (2) 

(ij) 

where cJ(q) is the creation (annihilation) operator at the lat- 
tice site i and t is the nearest-neighbor (NN) hopping am- 
plitude. The explicit single-particle Schrodinger equation 
derived from Eq. (2) is detailed in Appendix A. The en- 
ergy band structure of this system can be obtained from the 
Hamiltonian in momentum space H — £ k ^i/j^k, where 
= (cik,C2k,C3k,C4k,C5k) r - The non-zero components of 
the 5 x 5 matrix are given by (H^) TT f = exp(/k • v TT /), y/k o 
where v TT / is the vector connecting two NN sites T and t', 
which results in the Hamiltonian 
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(3) 

Diagonalizing H k yields the energy bands E(k) depicted in 
Fig. 1(b). 

Interestingly, the band structure displays two different 
regimes. At fillings / = 1/5 and / = 4/5, two indepen- 
dent Dirac cones located at K + = (§f> _ ^§^) an d K_ = 
( || , ) are present within the first Brillouin zone, which is 
the same as for the standard honeycomb lattice where Dirac- 
like dispersion relations effectively describe spin- 1/2 rela- 
tivistic fermions. At half-filling / = 1/2, a flat band is present 
at the tip of a single Dirac cone located at To = (0, 0). We find 
that the wavef unctions associated with the flat band have zero 
amplitude at the green sites (t = 1, 3, 5) illustrated in Fig. 1(a), 
which is compatible with the localization property expected 
from their infinite effective band mass 15 . This peculiar config- 
uration, involving a flat band and a single Dirac cone, is also 
present in the Lieb lattice and leads to an effective Hamilto- 
nian describing spin-1 relativistic fermions 10-12 . We see from 
Fig. 1 (b) that the ECH lattice indeed contains both spin-1 /2 
and spin-1 relativistic dispersion relations. These two singular 
and distinct regimes could be reached in a cold-atom realiza- 
tion by simply tuning the atomic filling factor. In order to 
further explore the distinction between these two relativistic 
regimes, we show the density of states p(E) in Fig. 1 (d). 
We find that around E w ±\^3t, p(E) behaves linearly which 
is expected for spin- 1/2 relativistic fermions. Around E « 0, 
the p (E) shows a linear behavior and a sharp peak as a conse- 
quence of the flat band 15 . As we will demonstrate below, the 
van Hove singularities located at E « ±y/2t constitute natu- 
ral boundaries separating the spin- 1/2 and spin-1 relativistic 



regimes. At this point, let us comment on the interesting fact 
that the additional (green) lattice sites, i.e., x — 1,3,5, do not 
destroy the relativistic properties stemming from the back- 
ground honeycomb lattice. They rather enrich the quantum 
properties of the lattice in a non-trivial manner by inducing 
new relativistic regimes at various fillings. 




Figure 1: (a) The ECH lattice is characterized by five sites per 
unit cell, denoted by % = 1, ... ,5. The three NN vectors of the un- 
derlying honeycomb structure (red sites) are given by ai = (1,0), 
a 2 = (-1/2, >/3/2), a 3 = (-1/2,-^/3/2) while the NN vectors of 
the ECH lattice by v M = a M /2, where (1 = 1,2,3. (b) The energy 
band structure E(k)/t of the ECH lattice hosts two different kinds of 
Dirac- Weyl fermions: spin- 1/2 at f= 1/5,4/5 and spin- 1 at / = 1 / 2. 
(c) The intensity profile obtained from the six-laser configu- 

ration in Eq. (1) that would create an ECH lattice for cold atoms, (d). 
The density of states (DOS) of the band structure illustrated in (b). 
The two relativistic regimes (spin- 1/2 and spin-1) are separated by 
van Hove singularities at E = ±y/2t indicated by the vertical dotted 
lines. 



III. THE SPIN-1/2 AND SPIN-1 RELATIVISTIC REGIMES 

To demonstrate the above assertion that the low-energy ex- 
citations around fillings / = 1/5,4/5 and / = 1/2, are in- 
deed Dirac- Weyl fermions of different kinds, we obtain the 
low-energy effective Hamiltonians describing these excita- 
tions around the band- touching points, i.e., K± at E = ±y/3t 
and ro at E = 0. To do so, we linearize H(k) near K± or To, 
and subsequently project onto the subspace associated with 
the two (for / = 1/5 and / = 4/5) or three (for / = 1/2) 
touching bands. This leads to 



h 



1/2 



= Vi/ 2 (p x ai+/yx2) 
= vi(p x Si+p y S 2 ) 



at 



at 



/= 1/5,4/5 
/=l/2 



(4) 
(5) 



where p = k - K±(Io) and v x / 2 = V3t/4, V\ = V?>t/2 are 
the Fermi velocities of the spin- 1/2 and spin-1 relativistic 
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fermions. We note here that Vi = 2v 1 / 2 is perfectly compati- 
ble with their spin-1 and spin- 1/2 nature. While around / = 
1/5,4/5, the Oi ? 2 are the usual Pauli matrices acting on the 
two-dimensional subspace associated with the two touching 
bands, the effective Dirac-Weyl Hamiltonian around / = 1/2 
features the 3 x 3 matrices Si^, which fulfill the angular mo- 
mentum commutation relation [Si,Sj] = iZijkSk- Such an effec- 
tive Hamiltonian has been shown to describe spin-1 massless 
Dirac-Weyl fermions, as recently discussed in Refs. 10-12 . The 
spin- 1/2 and spin-1 Dirac-Weyl fermions therefore do indeed 
coexist in the ECH lattice at different fillings. 

IV. THE "HONEYCOMB-KAGOME" DECOUPLING 

The ECH lattice has a bipartite structure which allows for 
an effective decoupling between the red (t = 2,4) and green 
(t = 1,3,5) sites illustrated in Fig. 1 (a). Note that while the 
red sites form the background honeycomb lattice, the green 
sites constitute a Kagome lattice. It turns out that the band 
structure depicted in Fig. 1 (b) can be deduced from the en- 
ergy spectra describing these two sublattices 20 . Indeed, as 
demonstrated in Appendix A, one finds that the five energy 
bands associated to the ECH lattice are directly obtained from 
the relations 

(E/t)(k) = ±y/e K (k)+2, (6) 

(E/t)(k) = ±y/e H (k) + 3, (7) 

where and 8k (k) are the energy bands related to the 

decoupled honeycomb and Kagome lattices, and where we 
assume that E ^ 0. The band structure 8k (k) is illustrated 
in Fig. 2 (a), which shows Dirac points at 8k = 1 and a flat 
band at 8k = — 2. Note that the dispersion relation is quadratic 
in the vicinity of the flat band. The band structure of the 
ECH lattice can then be entirely understood from the spec- 
trum 8k (k). From the relation (6), one obtains a flat band 
at E/t = -2 + 2 = and Dirac points at E/t = ±y/l+2. 
Furthermore, the quadratic dispersion of the Kagome lattice 
around 8k ~ — 2 leads to the conical intersection at E/t = 
(see Fig. 2 (b)). Therefore both the spin-1 and spin- 1/2 Dirac 
structures stem from the background Kagome lattice, formed 
by the green sites in Fig. 1 (a). However, it is worth em- 
phasizing that the Kagome lattice alone does not display a 
spin-1 Dirac structure, which highlights the richness of the 
bipartite ECH lattice with respect to its underlying honey- 
comb and Kagome lattices. Furthermore, we note that the 
spectrum associated to the honeycomb lattice £#(Jfc) does not 
contribute to the band structure E(k) in a significant man- 
ner. If we inde ed om it the flat b and at 8k = — 2, we find that 
±y/e K (k)+2 = ±^/8 H (k) + 3 (see Fig. 2 (c),(d)). 

The "honey comb-kagome" decoupling described by Eqs. 
(6) and (7) also explains the location of the van Hove singu- 
larities in Fig. 1 (d). The honeycomb lattice presents van 
Hove singularities 21 at £# = ±1, which lead to the four peaks 
?XE/t = ±7+1 + 3 and E/t = ±7-1 + 3 in Fig. 1 (d). The 
linear behavior of the DOS around E/t ~ ±y/3 and E/t ^ is 
also easily deduced from the conical intersections stemming 
from the Kagome lattice, as discussed above. 



Finally, we note that the ECH lattice has a bipartite nature, 
with Nk = 3 and Nh = 2 sites per unit cell, where Nk,h re- 
spectively denotes the number of green (i.e. T = 1,3,5) and 
red (i.e. T = 2,4) sites. Under such conditions, and since the 
tunneling only occurs between red and green sites, one can 
apply the theorem of Ref. 15 , which stipulates that a flat band 
necessarily exists in the energy spectrum and that its weight is 
given by Nk /5 — = 1 /5 in the normalized DOS. 



(a) 4 (b) 3 




Figure 2: Energy spectrum for the (a) Kagome and (c) Honeycomb 
lattices. Figures (b),(d) show the mapping (6)-(7), leading to the 
ECH lattice band structure. 



V. PERTURBATIONS 

We proceed by investigating the effects of various pertur- 
bations on the ECH lattice and the two different Dirac-Weyl 
fermions, highlighting their similar and different natures. We 
will in particular consider the effects of external magnetic 
fields, spin-orbit coupling and charge density waves. Such 
perturbations will not only provide interesting insights into the 
physical properties of the ECH lattice, but they will also offer 
the possibility to distinguish between the different relativistic 
species in an experiment. 



A. Synthetic magnetic fields 

When the ECH lattice is subject to a uniform magnetic flux 
<E> per plaquette, the energy bulk gaps form a fractal structure 
in the E — <E> plane, which consists of two Hofstadter-Rammal 
butterflies separated by a flat band 22 . The butterfly spectra 
of the honeycomb, Kagome, and lattices have been re- 
ported in Refs 20 ' 23 ' 24 . When the Fermi energy is exactly lo- 
cated inside such gaps, the Hall conductivity of the system 
is quantized. A typical sequence of Hall plateaus is shown 
in Fig. 3, for a reasonably small magnetic flux per plaquette 

~ 0.05. The Hall plateaus clearly evolve differently within 
the different regimes. Around / = 1/5,4/5, i.e. £p = ±\/3t 9 
the Hall plateaus feature the anomalous double step sequence 
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d^ 2 = ±2(N +l/2)e 2 /h, where N is an integer. In the vicin- 
ity of this spin- 1/2 regime, each Dirac fermion contributes to 
the Hall conductivity according to Obirac = ( e2 /h)/2, i.e., the 
system exhibits the so-called half-integer anomalous quan- 
tum Hall effect 25 ' 26 . This is also the case for spin- 1/2 Dirac 
fermions in graphene 27 . Around half-filling, i.e., = 0, 
one observes the characteristic sequence cr^ = ±NoNe 2 /h 
describing the quantum Hall plateaus for integer- spin Dirac- 
Weyl fermions 13 , where Afo is the number of Dirac points 
crossing the flat band. For integer spin Dirac-Weyl fermions, 
the absence of the half-integer anomaly leads to a zero Hall 
conductivity plateau. In Fig. 3, we see the characteristic 
zero Hall conductivity plateau for integer spin Dirac-Weyl 
fermions and Nd — 1 , which is in agreement with the fact that 
a single Dirac cone is present in this spin-1 regime. 

Interestingly, a sharp change of behavior occurs at £p = 
±y/2t located at the van Hove singularities present in the DOS 
(see Fig. 1 (d)), which constitute the boundaries between 

the <7j!/ 2 and cr^ sequences. We point out that in the stan- 
dard honeycomb lattice, the van Hove singularities constitute 
boundaries between relativistic and non-relativistic regimes 26 , 
which is very different from the result presented here. The 
Hall conductivity sequence obtained from an ECH lattice sub- 
ject to a uniform magnetic flux therefore combines the two 

Hall sequences a^ 2 and cr^ of spin- 1/2 and spin-1 Dirac- 
Weyl fermions respectively. Obtaining the Hall sequence 
Ch(£f), such as presented in Fig. 3, would provide a clear 
signature for the coexistence of spin- 1/2 and spin-1 Dirac- 
Weyl fermions in the ECH lattice. 

In a cold-atom framework, such a study would require the 
presence of a uniform synthetic magnetic field within the opti- 
cal ECH lattice. This difficult, but realistic, task would require 
to engineer Peierls phases exp(/0y) that accompany the hop- 
ping of the atoms along the links j, in such a way that the to- 
tal product of the phases along a plaquette yields fin — ^ 
where is the magnetic flux per plaquette. Such phases could 
be induced by means of Raman-assisted tunneling (see Ref. 28 
for a review of synthetic gauge fields for cold atoms), as re- 
cently demonstrated experimentally in Ref. 29 . Signatures re- 
lated to the Hall sequences could then be obtained from den- 
sity measurements, as discussed in Ref. 30 ' 31 . 

B. Spin-orbit coupling 

An intrinsic spin-orbit (SO) coupling term, 

H so = iho £ {dljXdtj)-a a pc] a Cjp, (8) 

has been introduced in Ref. 17 to predict the quantum spin Hall 
effect in a model described by a standard honeycomb lattice. 
Here Xso is the SO coupling strength, djj and d?- are two vec- 
tors connecting the NNN sites i and j, and a is the vector of 
Pauli matrices acting on the spin. This term opens a bulk gap 
in the energy spectrum that is associated to a non-trivial 
Z 2 index and hosts topologically protected helical edge states 
(i.e., counter-propagating edge states with opposite spin) 17 . 



This NNN SO term has been generalized to other lattices ex- 
hibiting spin- 1/2 relativistic fermions 16 ' 33 , and systematically 
leads to non-trivial Z 2 phases. In lattices featuring effective 
spin-1 fermions, such as the Lieb or ^3 lattices, the situation 
is more subtle. The NNN SO term leads to a trivial phase 
for the ^3 lattice but a non-trivial phase for the Lieb lattice 
9,10,34 Therefore, the effect of the SO term on the ECH lattice, 
in which spin-1 and spin- 1/2 excitations coexist, is a priori a 
non-trivial problem. Let us comment on the fact that the NNN 
hopping defined on the ECH lattice leads to a Kagome struc- 
ture. However, the path-dependent phases associated to the 
hopping in Eq. (8) generates a radically different spectrum, 
featuring non-trivial bulk gaps (cf. below). Finally, we men- 
tion that the SO term (8) could, in principle, be engineered in 
optical lattices (cf. Refs. 8 ' 35-37 ). 

First, we show how the SO term affects the low-energy the- 
ory describing the two kinds of Dirac-Weyl fermions. In this 
limit the effective Hamiltonians are 

h l k /2 = Vi / 2(k x Oi+kyO 2 )-—aXso03 (9) 
hi = Vi(k x S l +k y S 2 )-2V3aXsoS3 (10) 

where a = ± is the spin index. We note that Eq. (9) 
holds for all the spin- 1/2 Dirac species, namely for all K± 
at / = 1/5,4/5. Therefore, the SO term generates the same 
mass term for all the spin- 1/2 relativistic excitations, and thus 
opens bulk gaps at the four independent spin-1 /2 Dirac points. 
Two bulk gaps also appear in the vicinity of the spin-1 Dirac 
point (i.e., / = 1/2) with the flat band being preserved by 
the SO term. Our result indicates that the bulk gaps associ- 
ated with the spin-1 Dirac-Weyl fermion are much bigger than 
the gap associated with the spin 1/2 Dirac-Weyl fermions (cf. 
also Fig. 4). 

We now compute the four Z 2 indices v#, with N = 1 , 2, 3, 4 
associated with the four bulk gaps opened by Hso- Since the 
ECH lattice possesses inversion symmetry, the Z 2 topological 
invariant Vn associated with the Mh bulk gap can be easily 
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Figure 3: Hall conductivity Oh(^f) as a function of the Fermi energy 
for « 0.05. Vertical dotted lines indicate the location of van Hove 
singularities at E/t = ±y/2 (cf. Fig. 1 (d)). 
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Table I: Parity-eigenvalue pattern at the four ^-invariant momenta 
Fi for the four different occupied bands. All the bulk gaps are asso- 
ciated with a non-trivial Z 2 index v = 1 . 



evaluated through the formula 38,39 

rinwr»)=(-i) viv - (id 

i=0 m— 1 

In this expression, £ 2m (r t -) = ±1 is the parity eigenvalue 
associated with the 2m-th occupied energy band, which is 
evaluated at one of the four ^-invariant momenta k = T/. 
The latter can be expressed as r ; = qiW;/2 + q 2 ra;/2 with 
tz/ = {0,1} and m ; - = {0,1}, where qi = 2^/3(1,^^3) and 
q 2 = 27r/3(l, — a/3). Choosing the site T = 3 inside the unit 
cell as the center of inversion, the parity operator acts as 

^ 3 [yi (r) , y/2 (r) , y/3 (r) , 1//4 (r) , i/a 5 (r)] 

= [i/Ai(-r + ai-a 2 ),V4(-r + ai),i/A 3 (-r), 

VA 2 (-r-ai),i/A 5 (-r + a 3 -ai)], 

(12) 

where 1/^ (r) is the single-particle wavefunction defined at site 
T. The eigenstates of Ho+so(k = H), as well as the parity 
eigenvalues of the occupied bands, are determined numeri- 
cally, yielding the results presented in Table I. We find that 
the expression (11) gives (— l) v = —1 for each bulk gap, in- 
dicating that the Z 2 phases generated by the SO term are all 
non-trivial. Therefore, the spin-1 regime of the ECH lattice 
behaves similarly to the Lieb lattice 9 ' 10 ' 34 . 
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Figure 4: The energy spectrum E = E(k), for an ECH lattice subject 
to the Haldane-type SO coupling (8) with periodic boundary condi- 
tions along one spatial direction, shows the presence of helical edge 
states in the vicinity of both the spin- 1/2 and the spin-1 regime. 

To further confirm these results, we diagonalize the ECH 
lattice in the presence of the SO term and consider periodic 



boundary conditions along one spatial direction. In this cylin- 
drical geometry, the energy spectrum features helical edge 
states within the four bulk gaps predicted by the non-trivial 
Z 2 index 16,17,33 (see Fig. 4). We note that the dispersion re- 
lations E(k) associated with the edge states at half-filling are 
similar to those obtained from the Lieb lattice 10 ' 34 , further 
illustrating the similarity between the spin-1 regimes of the 
ECH and Lieb lattices. In the vicinity of the spin-1 /2 regime, 
i.e., / = 1/5,4/5, the dispersion relation of the edge states is 
similar to those obtained in the Kagome lattice 16 . These re- 
sults show that the SO term acts in a non- trivial way, both for 
the spin-1 and spin- 1/2 regimes, indicating that the ECH lat- 
tice presents striking similarities with both the Kagome and 
the Lieb lattices. However, the presence of a SO term does 
not allow us to distinguish between the spin-1 and spin- 1/2 
regimes, as they would both give rise to the same quantum 
spin Hall effect. Finally, we mention that edge states struc- 
tures, such as depicted in Fig. 4, could be probed in atomic 
systems through Bragg spectroscopy 40 ^ 2 . 



C. Charge density waves (CDW) 

In the honeycomb lattice, a staggered potential is known to 
open a trivial bulk gap at half-filling 17 . Such a perturbation 
acts as local chemical potentials \1a = — \Ib, which take oppo- 
site values at the sites A and B constituting the unit cell of the 
honeycomb. Such a perturbation has been generalized for lat- 
tices featuring J/ > 2 sites per unit cell, such as the Kagome 16 
and the decorated honeycomb lattice 33 . This charge-density- 
wave term is expressed as H^ DW = diag(/li, . . . \l ), which 
reduces to the honeycomb staggered potential for J/ = 2 and 
jUi = — jU 2 . For lattices with JV > 2 exhibiting spin- 1/2 Dirac 
fermions, it was shown that this CDW takes the form of an ax- 
ial gauge field in the effective low-energy Hamiltonian 16 ' 33 , 
namely a gauge potential A which has opposite sign at two 
independent Dirac points. 

Here we are interested in the fate of the spin- 1/2 and spin- 
1 Dirac points when a charge density wave term H^ DW = 
£diag(/li, . . . ,^5) is added in the ECH lattice, where the lo- 
cal chemical potentials /i T can be individually tuned. In the 
vicinity of the four spin- 1/2 Dirac points K±, we find that the 
low-energy terms corresponding to the CDW take the form 

*(pfi = -(A> 1 +A^ 2 -A^ 3 )+A /, f=^) 
Kv)cdw = (A>i+A^a 2 +A^ 3 )+A /, /=| (14) 

where A l x = - /isjZ/VS, A l y = -2ji 3 +jis)Z/12, A{ = 
(ju 2 - /i4)//4, A = (2jUi + 3jU 2 + 2/^3 + 3^ 4 + 2jU 5 )/12, and 
/ = ± refers to the two Dirac points K±. Therefore, when 
M2 = AU> i- e -> A^ = 0, and similar to the results reported for the 
Kagome lattice, we find that the CDW acts as an axial gauge 
field. In other words, the low-energy Hamiltonians 

%) 1/2 =I> 1/2 (£v-<Vv, (15) 

V 



where si* =A l v / v i/2 an d — ■> can be expressed in 
terms of a gauge potential A which has opposite sign at two 
independent Dirac points. In this case the effect of the gauge 
potential srf^- on the spin-1 /2 Dirac- Weyl fermions is to move 
the positions of their Dirac points inside the Brillouin zone. 
At a given filling, the displacement of the two Dirac points 
are opposite, since the gauge field is axial. Furthermore, these 
displacements are in opposite directions at fillings f = 1/5 
and / = 4/5 (cf. Eq. (14)). In other words, if two Dirac 
cones come closer at / = 1 /5 as the CDW is increased, the 
two cones at / = 4/5 will separate. When the CDW is suf- 
ficiently strong, the two approaching Dirac points at / = 1/5 
will annihilate each other 3 ' 31 ' 32 , while the two Dirac points at 
f = 4/5 will survive. This process allows the destruction of a 
pair of spin-1 /2 species at a given filling while preserving the 
others. 

In contrast, by expanding the Hamiltonian around the spin- 
1 Dirac point at To, we find that the low-energy form of the 
CDW perturbation is more involved than for the spin- 1/2 
regime, and that it cannot be simply interpreted as a gauge 
field. Indeed, the low-energy limit of the CDW term cannot 
be written as a superposition of the three angular-momentum 
matrices S x ^ z which do not form a complete basis for 3 x 3 
matrices. This interesting result indicates that the spin- 1/2 
and spin-1 regimes of the ECH lattice should react differently 
to the CDW, and thus providing a mechanism for distinguish- 
ing them in an experiment. 

The possibility to destroy and preserve the spin- 1/2 and 
spin-1 fermions individually, using the CDW perturbation, is 
appealing. Here we report a selection of relevant configura- 
tions that achieve this goal. 

• jii = M4 and Mi = M3 = M> In this case, the CDW 
acts as a staggered potential for the background honey- 
comb lattice. It destroys all the spin- 1/2 Dirac points 
at K± and / = 1/5,4/5 by opening trivial bulk gaps 17 . 
By setting \i\ = /I3 = /i 5 , the CDW do not perturb the 
localized states defined at the green sites (t = 1,3,5) 
of Fig. 1 (a). The flat band is therefore preserved. In 
addition, we find that when /I2 = — AU> the spin-1 Dirac 
point at To only survives for ^1,3,5 = 0. This situation 
is illustrated in Fig. 5 (a). 

• jl2 = M4 a nd ji\ = /I3 = /I5. In this case, A l x =A l y =A[ = 
and Aq 7^ 0, thus the spin- 1/2 Dirac points are simply 
shifted in energy. Since jl\ = /I3 = /I5 , the flat band is 
preserved but the spin-1 Dirac point at To is generally 
destroyed. This situation is illustrated in Fig. 5 (b). 

• fi2 = Ha and arbitrary In this case, the CDW acts 
as a non-trivial axial gauge field and the spin- 1/2 Dirac 
points move inside the Brillouin zone in opposite direc- 
tions. Therefore, for small CDW, the spin- 1/2 fermions 
are all preserved (cf. Fig. 5 (c)). For larger CDW, two 
fermions generally annihilate each other at / = 1 /5 or 
f = 4/5 (the displacements being in opposite direction 
at these fillings), in which case only one spin- 1/2 regime 
survives (cf. Fig. 5 (d)). In addition, for arbitrary /ii,3,5 
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Figure 5: Energy bands E = E(k x ,k y ) for different configurations of 
the CDW: (a) \ii = — JU4 = 1 and Jii 3 5 = 0: the spin- 1/2 regimes are 
gapped while the spin-1 regime is preserved, (b) 112 = M4 = 1 an d 
Ml, 3, 5 =0.5: the spin-1 regime is gapped while the spin- 1/2 regimes 
and the flat band are preserved, (c) M2 = M4 = 0-5, Ml = 0.25, 
jU3 = —0.6 and ^5 =0.1: the spin-1 regime is gapped, the flat band 
is perturbed, and the robust Dirac points at / = 1/5 move closer (cf . 
dotted circles), (d) The same CDW configuration as in (c) but mul- 
tiplied by a factor 1.4 (i.e. \i 2 = M4 = 0.7, jii = 0.35, fi 3 = -0.84 
and /I5 = 0.14 ): the Dirac points at / = 1/5 annihilate each other 
(cf. dotted circles), and only the spin- 1/2 Dirac points at / = 4/5 
survive. 



the flat band and the spin-1 Dirac fermion are generally 
destroyed. 



Therefore, by selecting the appropriate configuration of 
the CDW perturbation, one is able to engineer a system 
whose band structure displays zero, one or two spin- 1/2 Dirac 
regimes, a flat band or not, a spin-1 Dirac regime or not. 



VI. CONCLUSIONS 



We have investigated the ECH lattice which features both 
spin- 1/2 and spin-1 Dirac -Weyl fermions at different filling 
fractions. By using an intriguing mapping, we have shown 
that the underlying Kagome and honeycomb structures of the 
ECH lattice play a crucial role in determining the band struc- 
ture of the ECH lattice. We have also explored several types 
of perturbations on the lattice which offer a powerful control 
over this rich system. It is certainly tempting to extend this 
scenario to include collisional interactions between the spins. 
This will not only allow for exotic new phases 43,44 and appli- 
cations of models such as the Kitaev anyonic model 45 , but will 
also hopefully shed light on open questions at the forefront of 
condensed-matter physics. 
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Appendix A: The single-particle Schrodinger equation and the 
"honeycomb-kagome" mapping 

The single-particle Schrodinger equation describing non- 
interacting particles on the ECH lattice can be directly de- 



where ai = (1,0), a 2 = (-1/2, \/3/2), a 3 = (-1/2, -y/3/2) 
(see main text). For E ^ 0, one can decouple (Al) into two 



((E/t) 2 
((E/t) 2 
((E/t) 2 

and 



Writing \jf T (x) = exp(z* • x)(j) T , one finds the two separate 
eigensystems 



(<Pi 
e K (k) <fe 

and 



where we introduced the dimensionless quantities 
e K (k) = (E(k)/t) 2 -2, e H (k) = (E(k)/t) 2 -3. (A7) 



rived from the second-quantized Hamiltonian (2). Denoting 
the wave function at lattice site T = 1 , . . . , 5 by (x) , with 
spatial coordinate x, one finds the set of coupled equations 



(Al) 



subsets of equations describing the red (i.e., T = 2,4) and 
green sites (i.e., T = 1,3,5) independently. We find 



I 

The two decoupled systems, described by the Hamiltonians 
Hjj(k) and ##(*), correspond to the honeycomb and Kagome 



(E/t)\jfi(x) = y/4 (x + a 2 /2) + y 2 (x-a 2 /2), 

(E/t) y/ 2 (x - a 2 /2) y/i (x) + \y 3 (x-a 2 /2 + ai /2) + (* - a 2 /2 + a 3 /2) , 
(E/t) y/3 (x - a 2 /2 + a x /2) = y/ 4 (x - a 2 /2 + a x ) + \j/ 2 (x - a 2 /2) , 

(E/t)\i/ 4 (x-a 2 /2 + ai) = \j/i(x- a 2 + ai) + \j/ 3 (x- a 2 /2 + ai/2) + \j/ 5 (x-a 2 /2-a 3 /2 + ai), 
(E/t)\f/ 5 (x-a 2 /2-a 3 /2-\-a 1 ) = \\f 2 (x-a 2 /2 + a\ - a 3 ) + \j/ 4 (x - a 2 /2 + ai), 

I 



I 

- 2) yi (x) = y 3 (x + a 2 /2 - ai / 2) + i/a 3 (x - a 2 /2 + a : / 2) + 1// 5 (* + a 2 /2 - a 3 /2) + 1// 5 (* - a 2 /2 + a 3 /2) , (A2) 
-2)y/3(x-a 2 /2 + ai/2) yiW + + -*2) + ^(x-a 2 /2^a 3 /2) + i/A 5 (x-a 2 /2-a 3 /2 + ai), 
- 2) a 2 /2-a 3 /2 + a 1 ) = {x + a\ -a 2 ) + \j/\ (x + a\ - a 3 ) + \j/ 3 (x- a 2 /2 + a\ /2) + i// 3 (jc- a 2 /2 + 3a\ /2 - a 3 ), 

(A3) 

((E/t) 2 -3)\j/ 2 (x-a 2 /2) = \j/ 4 (x + a 2 /2) + \i/ 4 (x-a 2 /2 + ai) + \i/ 4 (x-a 2 /2 + a 3 ), (A4) 
((E/t) 2 -3)\\f 4 (x-a 2 /2 + a x ) = \jf 2 (x-a 2 /2) + Y 2 (x-3a 2 /2 + a 1 ) + Y 2 (x-a 2 /2 + ai-a 3 ). 

I 



cos*- (a 2 —a{)/2 cos*- (a 2 — a 3 )/2\ 1 0i 

: 2 ( cos*- (a 2 — ai)/2 cos*- (a 3 — fli)/2 I 3 

1 cos* - (fl 2 — a 3 )/2 cos*- (fl3 — fli)/2 / \ ^5 




(A5) 



£//(*) 





e -ika v q 



£)-*<*'(£ 



(A6) 
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lattices formed by the red (i.e., T = 2,4) and green (i.e., 
T = 1,3,5) sites, respectively. The energy bands of the two 
subsystems are given by 

e K (k) = (1 ± v^Ak-3), 8 K (k) = -2, (A8) 

£H = ±\t eihav l (A9) 
v=l 

where A k = cos 2 [k- (a2 — a3)/2] + cos 2 [k • (a3 - ai)/2] + 
cos 2 [k • (ai - a 2 )/2]. We note that ^/4A k - 3 = | J* =1 e* fl v j, 
thus we find that ±<y/e K (k)+2 = ±^£#(4) + 3 (for £ ^ 0). 



The latter result shows that the band structure describing the 
ECH lattice, E(k) ^ 0, can be equally obtained from the 
Kagome spectrum 8k (k) through the relations (A7)-(A8), or 
from the honeycomb spectrum through the relations 

(A7)-(A9). Besides, we note that the flat band at E = can 
also be deduced from the Kagome subsystem through (A7)- 
(A8), although we stress that this mapping is only strictly 
valid for E ^ 0. The band structures 8k ,//(£) are illustrated 
in Figs. 2 (a),(c). The five energy bands associated to the 
ECH lattice, obtained through the relations E/t = ±\/8k + 2 
and E/t = + 3, are depicted in Figs. 2 (b),(d). 
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